Dynamic Model of Super- Arrhenius Relaxation Rates in Glassy Materials 
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Super-Arrhenius relaxation rates in glassy materials can be associated with thermally activated 
rearrangements of increasing numbers of molecules at decreasing temperatures. We explore a model 
of such a mechanism in which string-like fluctuations in the neighborhood of shear transformation 
zones provide routes along which rearrangements can propagate, and the entropy associated with 
critically long strings allows the rearrangement to be distributed stably in the surrounding material. 
We further postulate that, at low enough temperatures, these fluctuations are localized on the 
interfaces between frustration-limited domains, and in this way obtain a modified Vogel-Fulcher 
formula for the relaxation rate. 



Recent developments in the shear-transformation-zone 
(STZ) theory of amorphous plasticity Q, Q{, especially 
its success in accounting for the behavior of bulk metallic 
glasses have prompted us to take a fresh look at the 

super- Arrhenius rates that characterize relaxation mech- 
anisms in such materials. In our STZ calculations so far, 
we have simply deduced these transition rates from mea- 
sured linear viscosities and used them in predicting, for 
example, nonlinear plastic responses to driving forces. To 
make further progress, it will be useful to have a deeper 
understanding of the physical mechanism that underlies 
these rates. (See Q for a recent summary of a wide va- 
riety of research in this and related fields.) 

The main assumption of the STZ theory is that a de- 
formable glassy material is fundamentally an elastic solid 
in which irreversible molecular rearrangements occur at 
special sites - so called "flow defects" or "STZ's." 0,00 
The STZ theory therefore differs from mode-coupling 
theories [HE! whose starting point is a liquid-like model. 
Our version of STZ theory differs also from the earlier 
flow-defect theories in that it models the defects as two- 
state systems that carry information about the internal 
state of the material. In this way, the theory describes 
both the onset of plastic flow at a yield stress and a va- 
riety of memory effects. 

In our most recent version of the STZ theory 0, we 
characterized the configurational de gre es of freedom, i.e. 
the inherent states of the system by an effective 

disorder temperature which, under noncquilibrium con- 
ditions, can be different from the temperature of the ther- 
mal bath. We supposed that the STZ's are especially de- 
formable local density fluctuations that are far out in the 
wings of the effective thermal distribution. Our implicit 
assumption was that only a narrow range of these config- 
urational fluctuations is sufficiently deformable yet popu- 
lous enough to participate in plastic deformation. We are 
concerned here with several different but closely related 
transformation rates pertaining to the STZ's. Specifi- 
cally, we want to compute the rate at which STZ's switch 
from one orientation to another during shear deforma- 
tions, the rates at which they are created and annihilated, 
and the rate at which the effective temperature relaxes 



toward the temperature of the heat bath. More gener- 
ally, we are interested in the mechanisms by which an 
amorphous solid makes transitions between its inherent 
states. 

As has long been recognized in theories of glass dy- 
namics (e.g. see Adam and Gibbs [HI), a qualitative 
reason why transition rates become anomalously slow 
at low temperatures is that, as the temperature de- 
creases, the statistically most probable transition states 
are those that involve increasingly large numbers of 
molecules. There have been several recent attempts to 
make this idea more precise. Notably, Xia and Wolynes 
[l^j l have postulated that the transitions are nucleated 
by liquid-like droplets characteristic of an incipient ran- 
dom first-order transition between solid-like and liquid- 
like (jammed and unjammed) phases. This postulate, 
plus some scaling analysis, produces the Vogel-Fulcher 
result. An alternative picture has been proposed by Gar- 
rahan and Chandler whose transitions are enabled 
by fluctuating patterns of mobile defects in an other- 
wise jammed system. The mobility of these defects - 
as opposed to our immobile STZ's - produces distinc- 
tive behavior in this model. In particular, Garrahan and 
Chandler conclude that relaxation times should obey the 
Arrhenius formula only at low temperatures rather than, 
as is more commonly supposed, changing from Arrhe- 
nius to super- Arrhenius as the temperature decreases. So 
far as we can tell, neither of these approaches attributes 
enough structure to the transition mechanism to describe 
the stress-driven change from thermally assisted viscous 
creep to superplastic flow that emerges from the present 
STZ theory. Our purpose here is to examine a transition 
mechanism that relates directly to our dynamical model 
of the STZ's. 

The model that we shall explore is motivated in part by 
work of Glotzer and colleagues 0, 0] , who discovered 
in molecular-dynamics simulations that transitions be- 
tween inherent states in glass-forming liquids take place 
via motions of string-like groups of molecules. We pos- 
tulate that, at temperatures low enough that most of the 
system is tightly jammed, localized molecular rearrange- 
ments might be cntropically enabled by strings of small 
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molecular displacements that distribute the disturbance 
throughout larger parts of the material. In granular ma- 
terials, our hypothetical mechanism might be visualized 
as a kinetic fluctuation that allows molecules to undergo 
small displacements along a force chain. 

For simplicity, consider first just the spontaneous STZ 
creation rate, that is, the rate at which STZ's are created 
by thermal fluctuations in the absence of external driv- 
ing. It is useful to think of STZ creation as an event in 
which the glassy analog of a vacancy and an interstitial 
first form, then move away from one another, and finally 
stabilize at an indefinitely large separation. More gener- 
ally, the formation of an STZ is a spontaneous increase 
in the configurational disorder of the system, as mea- 
sured by the intensity of density fluctuations. Suppose 
that, with a probability that we must calculate, the ma- 
terial in the neighborhood of this event contains a string 
of relatively loose molecules that provides a route along 
which the "vacancy" and the "interstitial" can propagate. 
Our strategy is to estimate the height of the free-energy 
barrier over which this system must fluctuate in order 
for it to become energetically favorable for the string to 
lengthen without bound. When that happens, we postu- 
late that a dynamically stable STZ has formed. In short, 
we propose to solve a nucleation problem where the reac- 
tion coordinate is the length of this string. The entropy 
associated with different string configurations is a mea- 
sure of the number of routes across this energy barrier 
and, therefore (see 0), reduces the free energy of the 
barrier for purposes of computing the nucleation rate. 

To describe the string model outlined in the preceding 
paragraphs, we let the string have length N in units of 
some characteristic molecular length, say, £, and suppose 
that it occupies a region of size R in the neighborhood 
of the emerging STZ. The total excess free energy of the 
system consists of several parts, which we denote: 

AG( A, R, T) = AGoo + N e - T S(N, R) + E mt (A, R) . 

(1) 

The first term, AGoo, is the bare activation energy for 
the transition, that is, the energy required to form the 
"vacancy" and the "interstitial." Until these two defects 
separate from each other, it will be energetically favor- 
able for them simply to recombine; thus, especially at low 
temperatures, we need the string to enable the system to 
reach a stable configuration. At high enough temper- 
atures, AGoo ought to become the ordinary Arrhcnius 
activation energy. 

The remaining terms on the right-hand side of Eq. Q 
describe the string, eo is the energy per unit step along it. 
S (A, R) is its entropy, which we obtain by computing the 
number of random walks of A steps extending a distance 
R. In the limit of large A, the number of such walks, say 
W(N,R), is approximately 



W(N, R) « constant x q N exp 



R 2 



2N£ 2 



(2) 



where q is the number of choices that the walk can make 
at each step. Thus, 



S(N,R) K.vk B N -k B 



R 2 



2N£ 2 



In q. (3) 



The last term in Eq.lJIJ, E int (N, R), is a repulsive in- 
teraction energy that accounts for the fact that no two 
pieces of the string can occupy the same position at the 
same time. This part of the analysis resembles Flory's 
calculation of excluded- volume effects in polymers. [18| 
Following Flory, we assume that Ei n t (A, R) is approxi- 
mately the square of the string density integrated over 
the volume occupied by the string. Therefore, using 
Flory's mean-field approximation, also in the limit of 
large A, we write 



Eint (N, R) w k B T int 



2 nd 



N 2 £' 



(4) 



where k B T n t is a repulsive energy (which contains di- 
mensionless geometric factors) and d is the dimensional- 
ity of the space in which this string exists. As we shall 
argue, it is not necessarily true that d = 3. 
Combining these terms, we have 



AG(N, R, T) AG C 



k f 



T 



R 2 



2N£ 2 



uN(T-T ) 

n 2 e d 

+ Tt, 



R d 



(5) 



where To = eo/{vk B ). The activation barrier is a saddle 
point in the A, R plane. That is, it is a minimum of 
AG(N, R, T) as a function of R (the smallest free energy 
for fixed A) and a maximum as a function of N (the high- 
est point along the reaction path). The two i?-dependent 
terms have a minimum at R = R* (JV, T) , where 



A :i 



[R*(N,T)}^k+—, 



(6) 



which is the Flory expression for the swelling of a d- 
dimensional polymer chain. Inserting this result into the 
independent terms in (J^J, we find that the activation 
energy has the following form as a function of A: 

AG*(N,T) = AG(N,R*,T) 

RJ AGoo + constant x T d /^ jv(4-«0/(<H*) 
- v N k B (T — To). (7) 

The second term on the right-hand side is positive and, 
for 1 < d < 4, is dominant for small enough A; the third 
term dominates at large A. For T > T , the activation 
energy goes through a maximum at N = A*(T), where 



N*(T) oc 



T d/(d+2) 



(T - To) 



(8) 



As in conventional nucleation theory, this fluctuation 
most probably will collapse for N < N*, but will grow 
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without bound if N becomes larger than N*. Thus the 
activation energy AG* (T) is the value of AG* (N, T) at 
its maximum, that is, 

AG*(T) = AG(N*,R*,T) 



T 2 ( d -v 

« AGoo + constant x 1=5— • (9) 

For the naively expected case of d = 3, these results 
are entirely unsatisfactory. The T-dependent factor in 
the activation energy, T 3 / 4 /(T - To) 1 / 4 , has too weak 
a divergence to be consistent with experimental data. 
Moreover, the energy scale is wrong. The implicit pic- 
ture is one in which the string consists of a chain of N 
monopolar, vacancy-like fluctuations, so that eo would 
be roughly equal to /i^ 3 , where /1 is the shear modulus 
and i is the molecular length scale introduced previously. 
Such an energy would be of the order of an electron volt, 
and would correspond to a temperature To in the range 
of 10 4 K - too large for our purposes by about two orders 
of magnitude. 

An apparently more plausible picture, and one which 
pertains specifically to the molecular structure of glassy 
materials, emerges from the concept of "frustration- 
limited domains," introduced by Kivelson et al. [l9ll20|. 
Their idea is that, in a glass-forming material, the en- 
ergetically preferred structure of small clusters of the 
constituent molecules is one that cannot tile an infinite 
space. That is, the energetically favorable short-range 
order is "frustrated" because it cannot extend over long 
distances. Thus a quenched glass may consist of many 
domains, inside of which the molecules have arranged 
themselves so as to have their preferred local coordina- 
tions - or some approximation thereto; but these coor- 
dinations are violated on the interfaces between the do- 
mains. Accordingly, we speculate that the STZ activity 
is localized on a network of two-dimensional interfaces 
that separate the domains. In addition to giving us a ra- 
tionale for choosing d = 2 in the preceding analysis, this 
hypothesis allows the energy eo to be much smaller than 
before, because the fluctuations are occurring in regions 
where the molecules already are more loosely bound to 
each other than they are within the bodies of the do- 
mains. 

Choosing d — 2 and restoring missing constants, we 
write Eqs. ©, ©, and @ as follows: 



[R*(N,T)} 2 



IT 



1/2 



l 2 A 3 / 2 ; 



A7*(T)«-L TmtT 



2v 2 (T-T ) 2 ' 



and 



AG*(T) AG C 



Tint T 



k B 2v(T-T )~ 



(10) 



(11) 



(12) 



This result exhibits the well-known Vogel-Fulcher linear 
divergence at T = To; therefore we know from earlier 
analyses (e.g. see [i^) that it will agree with experimen- 
tal data near the glass temperature. 

In the limit T — > 00, however, Ea. (|12fl predicts an 
excess, Arrhenius-like activation energy of magnitude 
k B Ti n t/2 v. (The situation is worse in three dimensions, 
where the activation energy grows like T 1 / 2 at high tem- 
peratures.) This physically unrealistic high-temperature 
behavior is a result of the fact that our large- N, mean- 
field approximations for the interaction energy and the 
entropy fail when N becomes small. In this limit, the 
string disappears and the interaction energy should van- 
ish accordingly; but our approximation says that the ra- 
tio N 2 /R 2 in Eq.(gJ goes to a constant. Note that the 
failure of the large-A approximation, by definition, oc- 
curs at the same temperature where the system switches 
from super- Arrhenius to simple Arrhcnius behavior. This 
transition region apparently is where the system also 
switches from being solid-like to liquid-like. 

The theoretical analysis in implies that the New- 
tonian viscosity t]n(T) can be written in the form 
ijoo exp [AG*(T)/fce T] up to slowly varying logarithmic 
corrections. In that analysis as well as here, the Arrhc- 
nius part of AG*(T), i.e. AGoo, is the STZ formation en- 
ergy; and the prefactor in the plastic strain rate, propor- 
tional to the Boltzmann factor exp (— AGoo/k B T), is the 
equilibrium density of STZ's. The non- Arrhenius part of 
AG*(T) in Ea. l|12l) determines the nonequilibrium rate 
factor denoted in [4j (up to a prefactor) by exp [— a(T)]; 
that is, a(T) » T int /\2u(T - T )]. Thus, in the spirit 
of Adam and Gibbs [12j, our super- Arrhenius rates are 
truly nonequilibrium quantities. They describe transi- 
tions between near-equilibrium, inherent states and not, 
as sometimes has been assumed, an equilibrium distribu- 
tion associated with the states themselves. For example, 
it is assumed in 0>E]>II2>IIii that the equilibrium density 
of flow defects has a Vogel-Fulcher form; and the Cohen- 
Grest model attributes super- Arrhenius behavior to 
percolation of liquid-like regions in equilibrated states. 
In our opinion, the nonequilibrium interpretation is the 
more natural of the two possibilities. 

In Fig.(l), we illustrate both the agreement near T 
and the asymptotic disagreement at higher tempera- 
tures by plotting AG*(T) obtained from our theory 
and from measurements of r]N(T) for the bulk metal- 
lic glass Zr4i.2 Tii3.s CU12.5 Niio Be22.5- We obtained the 
experimental points by first fitting the high-T (above 
T = 900 K) part of the data in jH E3| with the Ar- 
rhenius function t/n = rj^ exp (AGoo /&b T), finding 
r?oo = 1-78 x 10~ 10 and AGoo/k B T = 24300 K. We 
then followed the procedure described in |25j and plotted 
T log (77/7700) as a function of the temperature. To com- 
pare this data with the prediction of Ea. ((T2"|) . we used 
v = In 4, T = 515 K and T lnt = 2800 if. As expected, 
the theory fails above T = 700 K. It is interesting to 
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FIG. 1: Fits of metallic glass data from |26l l27 | . The dashed 
line is a fit of the low temperature data with our theory. Pa- 
rameters are AG x /k B = 24300 K, T = 515, T in t = 2800, 
and v — In 4. 

note, however, that the deformation measurements re- 
ported in [53 and discussed in 0,0] were made in the low- 
temperature regime, where Eq. jT^j) fits the data quite 
well. (Our theoretical comparisons to data for structural 
glasses such as those discussed in [2(| or [25| are qual- 
itatively the same as that shown here for the metallic 
glass.) 

Although the failure of Eq. (|12(1 at high temperatures is 
clearly related to a failure of our large- N approximations, 
we believe that simply improving the mathematics of our 
string model is not what is needed at this point. Rather, 
it seems to us that somehow we must construct a realistic 
model of the transition between solid-like and liquid-like 
glasses using physics that so far we have not brought to 
bear on this problem. The picture of a solid-like glass 
as a three dimensional mosaic of frustration-limited do- 
mains must break down at higher temperatures, where 
the domains must become smaller and the active, liquid- 
like regions between them must occupy a larger fraction 
of the system. As this happens, our strings - or whatever 
replaces them - must begin to look very different than 
they do in our simple two-dimensional approximation. 
At the most fundamental level, we believe we next must 
find a way to describe how a solid-like theory of the kind 
discussed here crosses over to a mode-coupling theory of 
a liquid-like glass forming material. 
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